Genomes in clinical care.

In the era of precision medicine, genome sequencing (GS) has become more affordable and the importance of genomics and multi-omics in clinical care is increasingly being recognized. However, how to scale and effectively implement GS on an institutional level remains a challenge for many. Here, we present Genome First and Ge-Med, two clinical implementation studies focused on identifying the key pillars and processes that are required to make routine GS and predictive genomics a reality in the clinical setting. We describe our experience and lessons learned for a variety of topics including test logistics, patient care processes, data reporting, and infrastructure. Our model of providing clinical care and comprehensive genomic analysis from a single source may be used by other centers with a similar structure to facilitate the implementation of omics-based personalized health concepts in medicine.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples.

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases.

PURPOSE: Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION). METHODS: PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation. RESULTS: A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants. CONCLUSION: GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Titin copy number variations associated with dominant inherited phenotypes.

BACKGROUND: Titinopathies are caused by mutations in the titin gene (TTN). Titin is the largest known human protein; its gene has the longest coding phase with 364 exons. Titinopathies are very complex neuromuscular pathologies due to the variable age of onset of symptoms, the great diversity of pathological and muscular impairment patterns (cardiac, skeletal muscle or mixed) and both autosomal dominant and recessive modes of transmission. Until now, only few CNVs in TTN have been reported without clear genotype-phenotype associations. METHODS: Our study includes eight families with dominant titinopathies. We performed next-generation sequencing or comparative genomic hybridisation array analyses and found CNVs in the TTN gene. We characterised these CNVs by RNA sequencing (RNAseq) analyses in six patients' muscles and performed genotype-phenotype inheritance association study by combining the clinical and biological data of these eight families. RESULTS: Seven deletion-type CNVs in the TTN gene were identified among these families. Genotype and RNAseq results showed that five deletions do not alter the reading frame and one is out-of-reading frame. The main phenotype identified was distal myopathy associated with contractures. The analysis of morphological, clinical and genetic data and imaging let us draw new genotype-phenotype associations of titinopathies. CONCLUSION: Identifying TTN CNVs will further increase diagnostic sensitivity in these complex neuromuscular pathologies. Our cohort of patients enabled us to identify new deletion-type CNVs in the TTN gene, with unexpected autosomal dominant transmission. This is valuable in establishing new genotype-phenotype associations of titinopathies, mainly distal myopathy in most of the patients.

UV-radiation and MC1R germline mutations are risk factors for the development of conventional and spitzoid melanomas in children and adolescents.

BACKGROUND: Genomic characterisation has led to an improved understanding of adult melanoma. However, the aetiology of melanoma in children is still unclear and identifying the correct diagnosis and therapeutic strategies remains challenging. METHODS: Exome sequencing of matched tumour-normal pairs from 26 paediatric patients was performed to study the mutational spectrum of melanomas. The cohort was grouped into different categories: spitzoid melanoma (SM), conventional melanoma (CM), and other melanomas (OT). FINDINGS: In all patients with CM (n = 10) germline variants associated with melanoma were found in low to moderate melanoma risk genes: in 8 patients MC1R variants, in 2 patients variants in MITF, PTEN and BRCA2. Somatic BRAF mutations were detected in 60% of CMs, homozygous deletions of CDKN2A in 20%, TERTp mutations in 30%. In the SM group (n = 12), 5 patients carried at least one MC1R variant; somatic BRAF mutations were detected in 8.3%, fusions in 25% of the cases. No SM showed a homozygous CDKN2A deletion nor a TERTp mutation. In 81.8% of the CM/SM cases the UV damage signatures SBS7 and/or DBS1 were detected. The patient with melanoma arising in giant congenital nevus (CNM) demonstrated the characteristic NRAS Q61K mutation. INTERPRETATION: UV-radiation and MC1R germline variants are risk factors in the development of conventional and spitzoid paediatric melanomas. Paediatric CMs share genomic similarities with adult CMs while the SMs differ genetically from the CM group. Consistent genetic characterization of all paediatric melanomas will potentially lead to better subtype differentiation, treatment, and prevention in the future. FUNDING: Found in Acknowledgement.

Clinical trio genome sequencing facilitates the interpretation of variants in cancer predisposition genes in paediatric tumour patients.

The prevalence of pathogenic and likely pathogenic (P/LP) variants in genes associated with cancer predisposition syndromes (CPS) is estimated to be 8-18% for paediatric cancer patients. In more than half of the carriers, the family history is unsuspicious for CPS. Therefore, broad genetic testing could identify germline predisposition in additional children with cancer resulting in important implications for themselves and their families. We thus evaluated clinical trio genome sequencing (TGS) in a cohort of 72 paediatric patients with solid cancers other than retinoblastoma or CNS-tumours. The most prevalent cancer types were sarcoma (n = 26), neuroblastoma (n = 15), and nephroblastoma (n = 10). Overall, P/LP variants in CPS genes were identified in 18.1% of patients (13/72) and P/LP variants in autosomal-dominant CPS genes in 9.7% (7/72). Genetic evaluation would have been recommended for the majority of patients with P/LP variants according to the Jongmans criteria. Four patients (5.6%, 4/72) carried P/LP variants in autosomal-dominant genes known to be associated with their tumour type. With the immediate information on variant inheritance, TGS facilitated the identification of a de novo P/LP in NF1, a gonadosomatic mosaic in WT1 and two pathogenic variants in one patient (DICER1 and PALB2). TGS allows a more detailed characterization of structural variants with base-pair resolution of breakpoints which can be relevant for the interpretation of copy number variants. Altogether, TGS allows comprehensive identification of children with a CPS and supports the individualised clinical management of index patients and high-risk relatives.

Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer.

BACKGROUND: Germline CDH1 pathogenic or likely pathogenic variants cause hereditary diffuse gastric cancer (HDGC). Once a genetic cause is identified, stomachs' and breasts' surveillance and/or prophylactic surgery is offered to asymptomatic CDH1 carriers, which is life-saving. Herein, we characterized an inherited mechanism responsible for extremely early-onset gastric cancer and atypical HDGC high penetrance. METHODS: Whole-exome sequencing (WES) re-analysis was performed in an unsolved HDGC family. Accessible chromatin and CDH1 promoter interactors were evaluated in normal stomach by ATAC-seq and 4C-seq, and functional analysis was performed using CRISPR-Cas9, RNA-seq and pathway analysis. RESULTS: We identified a germline heterozygous 23 Kb CDH1-TANGO6 deletion in a family with eight diffuse gastric cancers, six before age 30. Atypical HDGC high penetrance and young cancer-onset argued towards a role for the deleted region downstream of CDH1, which we proved to present accessible chromatin, and CDH1 promoter interactors in normal stomach. CRISPR-Cas9 edited cells mimicking the CDH1-TANGO6 deletion display the strongest CDH1 mRNA downregulation, more impacted adhesion-associated, type-I interferon immune-associated and oncogenic signalling pathways, compared to wild-type or CDH1-deleted cells. This finding solved an 18-year family odyssey and engaged carrier family members in a cancer prevention pathway of care. CONCLUSION: In this work, we demonstrated that regulatory elements lying down-stream of CDH1 are part of a chromatin network that control CDH1 expression and influence cell transcriptome and associated signalling pathways, likely explaining high disease penetrance and very young cancer-onset. This study highlights the importance of incorporating scientific-technological updates and clinical guidelines in routine diagnosis, given their impact in timely genetic diagnosis and disease prevention.

Detection and Characterization of a De Novo Alu Retrotransposition Event Causing NKX2-1-Related Disorder.

BACKGROUND: Heterozygous NKX2-1 loss-of-function variants cause combinations of hyperkinetic movement disorders (MDs, particularly childhood-onset chorea), pulmonary dysfunction, and hypothyroidism. Mobile element insertions (MEIs) are potential disease-causing structural variants whose detection in routine diagnostics remains challenging. OBJECTIVE: To establish the molecular diagnosis of two first-degree relatives with clinically suspected NKX2-1-related disorder who had negative NKX2-1 Sanger (SS), whole-exome (WES), and whole-genome (WGS) sequencing. METHODS: The proband's WES was analyzed for MEIs. A candidate MEI in NKX2-1 underwent optimized SS after plasmid cloning. Functional studies exploring NKX2-1 haploinsufficiency at RNA and protein levels were performed. RESULTS: A 347-bp AluYa5 insertion with a 65-bp poly-A tail followed by a 16-bp duplication of the pre-insertion wild-type sequence in exon 3 of NKX2-1 (ENST00000354822.7:c.556_557insAlu541_556dup) segregated with the disease phenotype. CONCLUSIONS: We identified a de novo exonic AluYa5 insertion causing NKX2-1-related disorder in SS/WES/WGS-negative cases, suggesting that MEI analysis of short-read sequencing data or targeted long-read sequencing could unmask the molecular diagnosis of unsolved MD cases. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Heterozygous UCHL1 loss-of-function variants cause a neurodegenerative disorder with spasticity, ataxia, neuropathy, and optic atrophy.

PURPOSE: Biallelic variants in UCHL1 have been associated with a progressive early-onset neurodegenerative disorder, autosomal recessive spastic paraplegia type 79. In this study, we investigated heterozygous UCHL1 variants on the basis of results from cohort-based burden analyses. METHODS: Gene-burden analyses were performed on exome and genome data of independent cohorts of patients with hereditary ataxia and spastic paraplegia from Germany and the United Kingdom in a total of 3169 patients and 33,141 controls. Clinical data of affected individuals and additional independent families were collected and evaluated. Patients' fibroblasts were used to perform mass spectrometry-based proteomics. RESULTS: UCHL1 was prioritized in both independent cohorts as a candidate gene for an autosomal dominant disorder. We identified a total of 34 cases from 18 unrelated families, carrying 13 heterozygous loss-of-function variants (15 families) and an inframe insertion (3 families). Affected individuals mainly presented with spasticity (24/31), ataxia (28/31), neuropathy (11/21), and optic atrophy (9/17). The mass spectrometry-based proteomics showed approximately 50% reduction of UCHL1 expression in patients' fibroblasts. CONCLUSION: Our bioinformatic analysis, in-depth clinical and genetic workup, and functional studies established haploinsufficiency of UCHL1 as a novel disease mechanism in spastic ataxia.

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank.

BACKGROUND: The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma. METHODS: We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping. RESULTS: Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11). CONCLUSIONS: At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.

Transcript-Specific Loss-of-Function Variants in VPS16 Are Enriched in Patients With Dystonia.

BACKGROUND AND OBJECTIVES: Our objective was to improve rare variant interpretation using statistical measures as well as publicly accessible annotation of expression levels and tissue specificity of different splice isoforms. We describe rare VPS16 variants observed in patients with dystonia and patients without dystonia, elaborate on our interpretation of VPS16 variants affecting different transcripts, and provide detailed clinical description of the movement disorder caused by VPS16 variants. METHODS: In-house exome and genome data sets (n = 11,539) were screened for rare heterozygous missense and putative loss-of-function (pLoF) variants in VPS16. Using pext (proportion expressed across transcripts) values from the Genome Aggregation Database (gnomAD), we differentiated variants affecting weakly and highly expressed exons/transcripts and applied statistical measures to systematically identify disease-associated genetic variation among patients with dystonia (n = 280). RESULTS: Six different heterozygous pLoFs in VPS16 transcripts were identified in 13 individuals. Three of these pLoFs occurred in 9 individuals with different phenotypes, and 3 pLoFs were identified in 4 unrelated individuals with early-onset dystonia. Although pLoFs were enriched in the dystonia cohort (n = 280; p = 2.04 × 10-4; 4/280 cases vs 9/11,259 controls; Fisher exact test), it was not exome-wide significant. According to the pext values in gnomAD, all 3 pLoFs observed in the patients with dystonia were located in the highly expressed canonical transcript ENST00000380445.3, whereas 2 of 3 pLoFs detected in 8 individuals without dystonia were located in the first exon of the noncanonical transcript ENST00000380443.3 that is weakly expressed across all tissues. Taking these biological implications into account, pLoFs involving the canonical transcript were exome-wide significantly enriched in patients with dystonia (p = 1.67 × 10-6; 4/280 cases vs 1/11,259 controls; Fisher exact test). All VPS16 patients showed mild progressive dystonia with writer's cramp as the presenting symptom between age 7 and 34 years (mean 20 years) that often progressed to generalized dystonia and was even accompanied by hyperkinetic movements and myoclonus in 1 patient. DISCUSSION: Our data provide strong evidence for VPS16 pLoFs to be implicated in dystonia and knowledge on exon resolution expression levels as well as statistical measures proved to be useful for variant interpretation.

Enhanced epigenetic profiling of classical human monocytes reveals a specific signature of healthy aging in the DNA methylome.

The impact of healthy aging on molecular programming of immune cells is poorly understood. Here, we report comprehensive characterization of healthy aging in human classical monocytes, with a focus on epigenomic, transcriptomic, and proteomic alterations, as well as the corresponding proteomic and metabolomic data for plasma, using healthy cohorts of 20 young and 20 older males (~27 and ~64 years old on average). For each individual, we performed eRRBS-based DNA methylation profiling, which allowed us to identify a set of age-associated differentially methylated regions (DMRs) - a novel, cell-type specific signature of aging in DNA methylome. Hypermethylation events were associated with H3K27me3 in the CpG islands near promoters of lowly-expressed genes, while hypomethylated DMRs were enriched in H3K4me1 marked regions and associated with age-related increase of expression of the corresponding genes, providing a link between DNA methylation and age-associated transcriptional changes in primary human cells.

Detection of mobile elements insertions for routine clinical diagnostics in targeted sequencing data.

BACKGROUND: Targeted sequencing approaches such as gene panel or exome sequencing have become standard of care for the diagnosis of rare and common genetic disease. The detection and interpretation of point mutations, small insertions and deletions, and even exon-level copy number variants are well established in clinical genetic testing. Other types of genetic variation such as mobile elements insertions (MEIs) are technically difficult to detect. In addition, their downstream clinical interpretation is more complex compared to point mutations due to a larger genomic footprint that can not only predict a clear loss of protein function but might disturb gene regulation and splicing even when located within the non-coding regions. As a consequence, the contribution of MEIs to disease and tumor development remains largely unexplored in routine diagnostics. METHODS: In this study, we investigated the occurrence of MEIs in 7,693 exome datasets from individuals with rare diseases and healthy relatives as well as 788 cancer patients analyzed by panel sequencing. RESULTS: We present several exemplary cases highlighting the diagnostic value of MEIs and propose a strategy for the detection, prioritization, and clinical interpretation of MEIs in routine clinical diagnostics. CONCLUSION: In this paper, we state that detection and interpretation of MEIs in clinical practice in targeted NGS data can be performed relatively easy despite the fact that MEIs very rarely occur in coding parts of the human genome. Large scale reanalysis of MEIs in existing cohorts may solve otherwise unsolvable cases.

New euchromatic variant dup(11)(p15.3p15.1) transmitted through two generations defined by low coverage whole genome sequencing.

We report on a 14-year old boy, his father, and his paternal uncle, all three carriers of a duplication of chromosomal region 11p15.3-p15.1. The aberration was transmitted by the grandmother, who is carrier of a balanced insertion 46,XX,ins(14;11)(q32.1;p15.3p15.1). In order to determine the precise molecular basis of this structural variant, we performed low-coverage whole genome sequencing on the boy's father. This approach allowed precise determination of the genomic breakpoints and revealed a duplication of 6.9 Mb, centromeric to the Beckwith-Wiedemann/Silver-Russell syndrome critical region in 11p15.5, that inserted in inverse orientation into 14q32.12 (according to HGVS nomenclature: NC_000014.8:g.92871000_92871001ins[NC_000011.9:g.12250642_19165928inv;T]). To our knowledge, this is the first report of a duplication of 11p15.3-p15.1 involving more than 40 genes and transmitted through two generations without apparent clinical effects.

Efficient and flexible Integration of variant characteristics in rare variant association studies using integrated nested Laplace approximation.

Rare variants are thought to play an important role in the etiology of complex diseases and may explain a significant fraction of the missing heritability in genetic disease studies. Next-generation sequencing facilitates the association of rare variants in coding or regulatory regions with complex diseases in large cohorts at genome-wide scale. However, rare variant association studies (RVAS) still lack power when cohorts are small to medium-sized and if genetic variation explains a small fraction of phenotypic variance. Here we present a novel Bayesian rare variant Association Test using Integrated Nested Laplace Approximation (BATI). Unlike existing RVAS tests, BATI allows integration of individual or variant-specific features as covariates, while efficiently performing inference based on full model estimation. We demonstrate that BATI outperforms established RVAS methods on realistic, semi-synthetic whole-exome sequencing cohorts, especially when using meaningful biological context, such as functional annotation. We show that BATI achieves power above 70% in scenarios in which competing tests fail to identify risk genes, e.g. when risk variants in sum explain less than 0.5% of phenotypic variance. We have integrated BATI, together with five existing RVAS tests in the 'Rare Variant Genome Wide Association Study' (rvGWAS) framework for data analyzed by whole-exome or whole genome sequencing. rvGWAS supports rare variant association for genes or any other biological unit such as promoters, while allowing the analysis of essential functionalities like quality control or filtering. Applying rvGWAS to a Chronic Lymphocytic Leukemia study we identified eight candidate predisposition genes, including EHMT2 and COPS7A.

Distinct Mutation Patterns Reveal Melanoma Subtypes and Influence Immunotherapy Response in Advanced Melanoma Patients.

The detection of somatic driver mutations by next-generation sequencing (NGS) is becoming increasingly important in the care of advanced melanoma patients. In our study, we evaluated the NGS results of 82 melanoma patients from clinical routine in 2017. Besides determining the tumor mutational burden (TMB) and annotation of all genetic driver alterations, we investigated their potential as a predictor for resistance to immune checkpoint inhibitors (ICI) and as a distinguishing feature between melanoma subtypes. Melanomas of unknown primary had a similar mutation pattern and TMB to cutaneous melanoma, which hints at its cutaneous origin. Besides the typical hotspot mutation in BRAF and NRAS, we frequently observed CDKN2A deletions. Acral and mucosal melanomas were dominated by CNV alterations affecting PDGFRA, KIT, CDK4, RICTOR, CCND2 and CHEK2. Uveal melanoma often had somatic SNVs in GNA11/Q and amplification of MYC in all cases. A significantly higher incidence of BRAF V600 mutations and EGFR amplifications, PTEN and TP53 deletions was found in patients with disease progression while on ICI. Thus, NGS might help to characterize melanoma subtypes more precisely and to identify possible resistance mechanisms to ICI therapy. Nevertheless, NGS based studies, including larger cohorts, are needed to support potential genetic ICI resistance mechanisms.

MDM2, MDM4 and EGFR Amplifications and Hyperprogression in Metastatic Acral and Mucosal Melanoma.

BACKGROUND: Mucosal and acral melanoma respond worse to immune checkpoint inhibitors (ICI) than cutaneous melanoma. MDM2/4 as well as EGFR amplifications are supposed to be associated with hyperprogression on ICI in diverse cancers. We therefore investigated the response of metastatic acral and mucosal melanoma to ICI in regard to MDM2/4 or EGFR amplifications and melanoma type. METHODS: We conducted a query of our melanoma registry, looking for patients with metastatic acral or mucosal melanoma treated by ICI. Whole exome sequencing, FISH and immunohistochemistry on melanoma tissue could be performed on 45 of the total cohort of 51 patients. Data were correlated with patients` responses to ICI and survival. RESULTS: 22 out of 51 patients had hyperprogressive disease (an increase in tumor load of >50% at the first staging). Hyperprogression occurred more often in case of MDM2/4 or EGFR amplification or <1% PD-L1 positive tumor cells. Nevertheless, this association was not significant. Interestingly, the anorectal melanoma type and the presence of liver metastases were significantly associated with worse survival. CONCLUSIONS: So far, we found no reliable predictive marker for patients who develop hyperprogression on ICI, specifically with regard to MDM2/4 or EGFR amplifications. Nevertheless, patients with anorectal melanoma, liver metastases or melanoma with amplified MYC seem to have an increased risk of not benefitting from ICI.

Cancer immune control needs senescence induction by interferon-dependent cell cycle regulator pathways in tumours.

Immune checkpoint blockade (ICB)-based or natural cancer immune responses largely eliminate tumours. Yet, they require additional mechanisms to arrest those cancer cells that are not rejected. Cytokine-induced senescence (CIS) can stably arrest cancer cells, suggesting that interferon-dependent induction of senescence-inducing cell cycle regulators is needed to control those cancer cells that escape from killing. Here we report in two different cancers sensitive to T cell-mediated rejection, that deletion of the senescence-inducing cell cycle regulators p16Ink4a/p19Arf (Cdkn2a) or p21Cip1 (Cdkn1a) in the tumour cells abrogates both the natural and the ICB-induced cancer immune control. Also in humans, melanoma metastases that progressed rapidly during ICB have losses of senescence-inducing genes and amplifications of senescence inhibitors. Metastatic cells also resist CIS. Such genetic and functional alterations are infrequent in metastatic melanomas regressing during ICB. Thus, activation of tumour-intrinsic, senescence-inducing cell cycle regulators is required to stably arrest cancer cells that escape from eradication.